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We analyze the charge transport through a single level quantum dot coupled to a normal (N) and 
superconducting (S) leads where the electron pairs exist either as the coherent (for temperatures 
below T c ) or incoherent objects (in a region T c < T < T*). This situation can be achieved in practice 
if one uses the high T c superconducting material where various precursor effects have been observed 
upon approaching T c from above. Without restricting to any particular microscopic mechanism we 
investigate some qualitative changes of the nonequilibrium charge current caused by the electron 
pair coherence. 



I. INTRODUCTION 

It is a firm experimental fact that the phase transi- 
tion from normal to superconducting states of the under- 
doped high T c materials is accompanied by appearance of 
a pseudogap .1} . While approaching the critical tempera- 
ture T c from above the single particle states are gradually 
depleted in a certain energy region \lu \ < A pg around the 
Fermi level. This phenomenon is often interpreted theo- 
retically as a precursor of the true superconducting gap 
which is usually predicted by the BCS-like treatments for 
temperatures T < T c . On a microscopic basis the pseu- 
dogap can be assigned to presence of the electron pairs. 
Above T c their long-range coherence is destroyed by the 
strong quantum fluctuations arising either from the re- 
duced dimensionality, or due to close neighborhood of 
the Mott insulating state or because of a competition 
with various kinds of ordering in the system. The inco- 
herent electron pairs have been unambiguously detected 
experimentally above T c in measurements of the large 
Nernst coefficient [2] and in observation of the Berezinski- 
Kosterlitz-Thouless phase fluctuations [3| . There is how- 
ever a great amount of controversy regarding the temper- 
ature extent where the incoherent pairs eventually exist. 
According to available experimental data their existence 
is established for at least a dozen Kelvin above T c but 
such region can perhaps spread over a much wider regime 
up to T* at which the pseudogap finally closes. 

Various tunneling techniques have been used for a long 
time as a useful tool for probing the single particle spectra 
of the correlated systems. Recent technologicalprogress 
of the spectroscopic methods such as the STM [4| , the k- 
resolved photoemission (ARPES) [5], the Andreev-type 
techniques [6] and the Fourier transformed scanning tun- 
neling spectroscopy [7( allow for a precise measurements 
of the energy, momentum and space dependent density 
of states. In the present work we propose to consider the 
N-S junction with the quantum dot located in between. 
In this setup one would be able not only to detect the 
(pseudo)gap in the single particle excitation spectrum 
but, due to activation of additional transport channels 
(mainly the Andreev reflections), there would be a pos- 
sibility to analyze the electron pair coherence. 



Our discussion here is not limited to any particular 
microscopic model describing the electron pair forma- 
tion and their eventual coherence. On rather general 
grounds we investigate the proximity effect which leads to 
a particle- hole mixing at small energies < A and fur- 
thermore we analyze its influence on the effective charge 
transport through the quantum dot. In section II we 
briefly introduce the problem, in particular explaining 
how we treat the coherent and incoherent electron pairs. 
Next, we show the QD spectrum without any correlations 
(section III) and for the limit of very strong correlations 
(section IV). The main part of our study is in the sec- 
tion V, where we discuss the differential conductance as 
a function of the applied bias V and temperature T, both 
below and above T c . 



II. FORMULATION OF THE PROBLEM 

For a description of the quantum dot (QD) coupled to 
a normal (N) and superconducting (S) leads we consider 
the single impurity Anderson model 



H = H 



N 



rf lck/3 CT + h.i 



(1) 



k,<r /3={N,S} 



Operators d a (d£ ) annihilate (create) the QD electrons 
whose energy is £d- The Coulomb potential U > de- 
scribes a repulsion between electrons of opposite spin a 
and the hybridization V£g is responsible for transferring 
electrons from the QD the normal (/? = N) or supercon- 
ducting (S) leads. 

We assume that the normal electrode is described 
by the Hamiltonian of noninteracting fermions Hn = 

Sk, CT {EltN-V'N)q IUTN CkoN- For description of the su- 
perconducting lead we use a general expression 
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where the two-body term V pa iring induces either the co- 
herent (for T <T C ) or incoherent electron pairs (above 
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T <T*. In such a case the pseudogap affects the single 
particle spectrum only through the BCS-like shape of the 
diagonal part Q 



G5(k,w) 



(7) 



FIG. 1: Temperature dependence of the single particle gap 
A(T) (solid line), the damping rate 7(T) (dashed line) and 
the superconducting order parameter A SC (T) (dotted line). 



T c ). 

Without specifying V pa iring nor restricting to any par- 
ticular microscopic mechanism of superconductivity we 
use the basic BCS-type results. For the superconducting 
state (below the critical temperature T c ) we introduce, in 
a standard way, the following retarded Green's function 
expressed in the Nambu representation 



GS(k,w)- 



Ak 



+ 6cS 



(3) 



where £ k/ 3 = £kfj ~ Hp measures the energies from the 
chemical potential. The off-diagonal terms Ak have as 
usually a meaning of the gap in the single particle exci- 
tation spectrum of S electrons. 

Since the excitation gap is known to develop well above 
the transition temperature we follow the arguments of 
Levin et al [8j and impose the following phenomenological 
Ansatz 



A k = A k . sc 



A k, P3 - 



(4) 



The first part in (g)) is related to the superconducting 
order parameter (c_kj,sCk|s) while the latter one comes 
from the pseudogap. For a quantitative discussion we 
propose the temperature dependence given by 



A k , sc (T) 



A k (0) 




for T < T c , 
for T >T C . 



(5) 



The pseudogap contribution A k , pg to the effective gap 
origins from the preformed electron pairs. Above T c 
their long range coherence is missing (hence a name of the 
incoherent pairs) and such pairs ultimately dissociate at 
a certain temperature T* . In analogy to ([5]) we propose 
to consider the following temperature dependence 



A k (T) = A k (0) Wi- 



fe) 



Absence of the off-diagonal long range order above T c 
implies that the retarded Green's function §3§ must re- 
duce to a diagonal structure for all temperatures T c < 
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with the quasiparticle dispersion = 

the usual coherence factors u^. = | (l -r -g— j — ± — ^ k 

This sort of behavior (J7J) can be derived on a microscopic 
level investigating the pairing interactions beyond the 
mean- field BCS framework Q. 

In addition to ([5]) and (O we have chosen for compu- 
tational purposes some phenomenological damping rate 
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0+ 
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for T < T c , 
for T r <T<T* 
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and assumed its momentum variation in the following 
way 7 k = 7 2 / (7 + -^j) where parameter 7 = 7 kF . Fig- 
ure (HJ) illustrates the temperature dependences described 
in this section. 



III. THE UNCORRELATED QD 

To introduce the formalism of our calculations we first 
start analyzing the equilibrium case \xl= Hr- In the stan- 
dard Nambu notation we express the retarded Green's 
function of the QD through the Dyson equation 



GSM 



(9) 



with two contributions to the matrix selfenergy YT^u) = 
Y, r N (u>) + Eg(w). This problem can be solved exactly in 
a case of the noninteracting QD (U = 0) when 



uj — Ed + i0 + 
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and the selfenergies simplify to [ic| 
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The BCS coefficients it k v k = assure that a proxim- 

ity effect in the QD occurs only for temperatures T <T C . 
Otherwise, the pseudogap appearing in the spectrum of S 
electrons affects the Green's function ([7]) via the diagonal 
parts of (Ql 
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FIG. 2: Spectral function pd{u) of the QD for (7 = with the 
right h.s. electrode being in a superconducting (solid line) and 
in a normal state (dashed line) . We used the isotropic energy 
gap A k = 0.173 and e d = -0.0573, T = 0.0573, 7 = 0.0173 
and set the half-bandwidth 73 as a unit for energies. 
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FIG. 3: The real and imaginary parts of the matrix selfenergy 
for the diagonal Ejj^w) (top panel) and off-diagonal terms 
£12 (w) (bottom panel). We used the same set of parameters 
as in figure [1] 




FIG. 4: The ground state spectral function /9d(w) of the QD 
for [7 = versus a varying position of the energy level Ed- One 
can note a clear particle-hole mixing (two Lorentians built 
around ±ed). 



trum. 

(i) Since S electrons can occupy no states in the en- 
ergy gap \lo\ < A thereby the line broadening gets 
reduced by 50% and in consequence the QD peak 
around Ed becomes more pronounced. 



(ii) A large amount of S electron states is cumulated 
near to = ±A (i.e. at the square root divergences 
in the density of the S lead). Efficiency of the 
hybridization VkS is there considerably enhanced 
leading to a depletion of the QD states at w = ±A. 



(iii) A role of well defined quasiparticles in the super- 
conducting state is played by the electron pairs. 
Due to the hybridization Vks such a particle hole 
mixing is also transfered onto the QD spectrum 
(the proximity effect). Indeed, in figures and 0] 
we notice that besides the Lorentzian peak centered 
around Ed there also appears its tiny mirror reflec- 
tion at — Ed- 



The hybridization coupling V^p can be conveniently 
replaced by the weighted density function Tp(oS) — 
27r^ k |Vk/3| 2 5(w— £\np)- For both electrodes being normal 
the QD spectral function Pd{u) = — ^G dll (u>) acquires 
a Lorentzian shape centered around the single particle 
level Ed (see the dashed line in figure [2]) with the effective 
broadening r = rAr+r s = 2rAr. 

If one of the electrodes is a superconductor with an 
isotropic (k independent) energy gap we can notice some 
qualitative (and quantitative) differences in the QD spec- 



To provide the arguments for the above mentioned ef- 
fects we plot in figure [3] the diagonal and off-diagonal 
parts of the matrix selfenergy £^(0;). An odd symme- 
try of Im [Ei2(w)] gives a non-vanishing Re [£3.2 (w)] f° r 
all energies located inside the energy gap \u>\ < A. In 
diagonal term En(u) the imaginary part is even (and 
negative) while the real part is odd (therefore vanishing 
at u> = 0). Similar quantitative behavior of the matrix 
selfener gy Ed(y ) has been previously reported by several 
authors [12l 0, [3, EH 1 however, no clear evidence of the 
particle-hole mixing has been emphasized so far. 
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FIG. 5: Spectral function Pd(tt;) of the QD in the limit U = oo 
obtained for Ed = — 0.0513, Tp = 0.01D assuming the isotropic 
energy gap A k (T = 0) = O.ID and T* = 2T C . 



IV. THE STRONG CORRELATION LIMIT 

It is well known both from the theoretical [ll| and 
experimental studies [H, [TtJ that the Coulomb interac- 
tions have a remarkable influence on transport properties 
through the QD. In particular, such correlations are re- 
sponsible for the Coulomb blockade (observed by oscilla- 
tions of the differential conductance) and, at sufficiently 
low temperatures, produce the Kondo resonance leading 
to enhancement of the conductance to the unitary limit 
value 2e 2 /h. 

In this section we consider the correlations focusing on 
the extreme limit of U = oo. Under such condition no 
double occupancy of the QD is allowed and one expects 
it to have a tremendous effect on the charge tunneling es- 
pecially in the anomalous channels involving the electron 
pairs. 

Excluding the doubly occupied states from the Hilbert 
space can be formally achieved using the auxiliary fields 



(13) 



where the boson 6^) anc i fermion operators f a corre- 
spond to an annihilation (creation) of the empty and 
singly occupied states on the QD. These new fields must 
obey the local constraint b'b + flf a = 1. 

There are various methods to deal with the local con- 
straint. For simplicity, we apply here the technique pro- 
posed by Le Guillou and Ragoucy [H[ where projecting 
out of the doubly occupied states is achieved by appro- 
priate commutation relations between the operators of 
auxiliary fields. For the present context ([T]) some neces- 
sary technical details have been previously discussed in 
Ref. [3. 

In the limit U = 00 the Dyson equation §§§ can 
be solved using the renormalized propagator g r d {u>) — 
(1 — n^ a ) g r ®(uj) and the matrix sclfenergy 

££H = [E°» + ££»] / (1 - n_„) (14) 

where the contribution YPp of noninteracting electrons is 
given in (fTTl [T3i|) . The other contribution Ei r originates 



from the correlations and under appropriate conditions 
leads to the Kondo effect [ll|. One finds [l(| 



(15) 



where T3 is the Pauli matrix and rik/3 denotes an average 
occupancy of the k-momentum in the /3-th lead given by 



1 



exp(|^ 



1 



for (3 = N, 
for (3 = S. 



(16) 



In figure [5] we show the spectral function pd(w) cal- 
culated for several temperatures in both, the supercon- 
ducting and pseudogap states. In a comparison to the 
previous situation U — we can notice that: 

(i) for temperatures T < T c there are visible two 
Lorentzian peaks, however their positions are more 
shifted from ±£d because of a finite real part of the 
matrix selfenergy HH), 

(ii) for very low temperatures (T < Tk) there appears 
a narrow Kondo resonance at the Fermi energy as- 
sociated with the spin singlet made of the QD and 
itinerant electrons (see the inset in Fig. [5]), 

(iii) in the pseudogap regime above T c we no longer 
observe a tiny Lorentzian at lo ~ — Ed and simul- 
taneously a dip of the spectral function at \u)\ = 
±A pg (T) gets smeared because of the damping ef- 
fects. 

For the particular set of parameters = 0.01D, Ed 
used in figure [5] we estimate the Kondo peak disap- 
pears for temperatures higher than Tk ~ 0.15T c .On the 
other hand, for temperature exceeding T* the QD spec- 
trum evolves back to its single Lorentzian peak centered 
around Ed- 



V. TRANSPORT PROPERTIES 

In order to study the noncquilibrium physics we use the 
Keldysh formalism. Applying a bias V leads to imbalance 
of the chemical potentials /ijv - = eV which induces the 
charge current J(V) = ~ej- t J2k A c kNa c *N<x) through 
the QD. 

Following the procedure described previously [Io| we 
express the charge current J(V) in terms of the following 
contributions 



J = Jll + J 12 + J22 + J A- 



(17) 



The first three components in equation (jTTJ) have the 
Landauer-type structure 

Jij(V) = y J da; Tij{V) [/(w - eV) - f(cj)} (18) 
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FIG. 6: The differential conductance G(V, T) versus the ap- 
plied bias V for a representative set of temperatures in the su- 
perconducting region (the lines) and for the pseudogap phase 
(the symbols). We used the same set of parameters as in fig- 
ure [5]and conductnce is expressed in units of Go = G(0, T*). 
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FIG. 8: The differential conductance Ga{V, T) of the Andreev 
current for a number of temperatures (see the legend). The 
zero bias enhancement (central peak) is due to the Kondo 
resonance. Above T c the Andreev current vanishes [f|. 




FIG. 7: Contributions of the normal Gn (the red line) and 
the anomalous channels G12 (blue), G22 (green), Ga (indigo) 
to the total differential conductance G(V,T). Temperature 
is T = 0.01T C (i.e. T < Tk) and other parameters have the 
same values as in figure [S] 

where the corresponding transmittance are defined by 

Tn(V) = -ImE^s |G n | 2 T N (u), (19) 
T 12 (V) = -2ImS^ 2)S Re [G U G* 12 ] T N (cu), (20) 

T 2 2{V) = -Im££ 2)S \G 12 \ 2 T n (lo). (21) 
The last contribution describes the Andreev current 

Ja(V) = \\^ t a(V) [f(u - eV) - f(uj + eV)] (22) 

with 

T A {V) = -ImE^ |G 12 | 2 T N (u). (23) 

This kind of current arises when electron from the N 
lead is converted into the Copper pair in the S elec- 
trode and simultaneously a hole is reflected back to the 
N lead. For a detailed discussion of such anomalous 
Andreev current see for instance the recent review ar- 
ticle [6J. In a case when the both leads are normal (i.e. 
A — * 0) there survives only Ju(V) current and more- 
over the transmittance lfl9]) simplifies to the usual form 



In figure [5] we plot the differential conductance 
G(V,T) = dJ{V)/dV as a function of the external bias 
V for several representative temperatures. In the super- 
conducting state we clearly notice a strong suppression 
of the charge current at small voltages \eV\ < A(T). 
Due to energy gap in the spectrum of S electrons the 
only possible process of a charge tunneling occurs then 
through the Andreev channel. At temperatures below T c 
the non- vanishing conductance G(V,T) ~ Ga(V, T) for 
\eV\ < A(T) is almost an order of magnitude smaller 
than the normal state conductance G(V,T*). 

For sufficiently low temperatures T <Tk 1 we observe 
formation of the Kondo peak (see figure E|) which af- 
fects the differential conductance at small bias \V\. In 
the present it gives some enhancement of the zero bias 
conductance (see the inset of figure [9]) through the An- 
dreev scattering. This zero bias anomaly is however 
rather residual as compared to the N-QD-N junctions 
[ill ]. A strong suppression of the zero bias ano maly 
has been previously pointed out by several authors [12j. 
We checked that the low temperature differential con- 
ductance G(0, T) fits very well the universal parabolic 
variation characteristic for the Kondo regime. 

For higher voltages, exceeding A(T), the dominant 
role to the charge transport comes from the normal cur- 
rent Jn(V). In figure [7] we show each of the contribu- 
tions to the total conductance for a small temperature 
T <T K <T C . The anomalous channels J i2 (V) and J22 (V) 
get activated outside the energy gap. But these contri- 
butions, as well as Ga(V), do quickly diminish for an 
increasing bias |V|. 

The in-gap conductance arising from the Andreev cur- 
rent is very sensitive to the temperature (see figure [5]). 
Already for temperatures higher then T ~ Tk the Kondo 
peak is washed out which gives a concomitant disappear- 
ance of the zero bias anomaly. Upon further increasing 
the temperature there occurs a gradual suppression of the 
Andreev current which completely vanishes at T — > T c 7 . 

For temperatures above T c the charge current is trans- 
mitted only via the normal Ju(V) channel. With an 
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FIG. 9: Temperature dependence of the zero bias differential 
conductance G(V = 0, T) for the set of parameters used in 
figure [5] Inset shows the enhancement at low temperatures 
which is due to the Kondo resonance. 

increasing temperature the pseudogap becomes gradu- 
ally filled in (see figure [6]) therefore the zero bias con- 
ductance smoothly increases, reaching a local maximum 
at T* . From basic considerations [llj it is well known 
that in the normal state (for T>T*) the differential con- 
ductance exponentially decreases with respect to T. We 
would like to stress that the non-monotonous tempera- 
ture variation of the zero bias differential conductance 
G(0,T) (shown in figure [9]) could be a useful method for 
identifying the temperature region of the incoherent elec- 
tron pairs. 

VI. CONCLUSIONS 

We have investigated the charge tunneling from the 
normal (N) to superconducting (S) leads via the quan- 
tum do t (QD) . In agreement with the previous studies 
[12, El, EH HH we find that the long range off-diagonal 
superconducting order induces the proximity effect in the 
quantum dot (independently of the on-dot Coulomb in- 
teraction U) . Besides the main Lorentzian peak near the 
single particle energy level Sd there appears another small 
counter-peak around —Ed due to the particle-hole mixing. 
Such physical situation occurs for all temperatures below 
T c . This effect could be practically observed in measure- 
ments of the differential conductance upon varying the 
gate voltage. 

In a presence of the strong on-dot correlations U there 
arise some additional qualitative features. One of them is 



the Kondo state which occurs at sufficiently low temper- 
atures, usually below IK [la]. Spectral function of the 
QD develops then a narrow peak at the Fermi energy. It 
has a qualitative influence on the low energy transport 
properties. Since the normal channel Jn(V) is forbidden 
for |eV| < A (because of the energy gap of S electrons) 
the only possibility for the charge transfer is through the 
Andreev tunneling. A magnitude of the Andreev cur- 
rent is by far smaller than the normal charge current but 
nevertheless, in analogy to the N-QD-N Kondo anomaly 
fill Efjj ] , we find some enhancement of the zero bias con- 
ductance for the N-QD-S setup. This anomaly shows up 
at temperatures T < Tk but the conductance does not 
reach the unitary limit value 2e 2 /h. 

We have also extended our study on the pseudogap 
state and found that incoherent electron pairs have a re- 
markable influence on the charge transport. Above T c the 
zero bias differential conductance G(0, T) systematically 
increases upon increasing temperature because the pseu- 
dogap is gradually filled in. Finally, at T*, it reaches a 
local maximum and for higher temperatures (in the nor- 
mal state) has a changeover to the exponential decrease 
versus T. Such nonmonotonus variation of the zero bias 
conductance (figure O can serve as a sensitive method 
to identify the temperature region T c < T < T* of the 
incoherent electron pairs. 

Some other effects which might eventually come into 
the play. One of them is a strong anisotropy of the gap 
function Ak well known for all the HTSC cuprates. To 
some extent this issue has been recently addressed by 
ETH group [19] . Authors investigated the conductance of 
charge current flowing between the conducting STM tip 
and the quasi-2D superconducting Cu02 planes involving 
the apical oxygen atoms (to be regarded as the QD). 
They focused on studying the inelastic scattering caused 
by the oxygen atoms' vibrations. In the future work we 
shall discuss the normal and anomalous currents for such 
anisotropic gap Ak of the d-wave symmetry. As known 
from the literature the gapless superconductivity does 
not preclude appearance of the Kondo resonance [H, HH . 
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